Critical quasiparticle theory: Scaling, thermodynamic and transport properties 
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We use the recently developed 0] critical quasiparticle theory to derive the scaling behavior 
associated with a quantum critical point (QCP) in a correlated metal. This is applied to the 
magnetic-field induced QCP observed in YbRli2Si2 (YRS), for which we also derive the critical 
behavior of the specific heat, resistivity, Griineisen coefficient, and the thermopower. The theory 
accounts very well for the experimental results. 



Introduction 

Recent advances in low-temperature experimental 
techniques have stimulated much interest in quantum 
critical phenomena, which comprise phase transitions at 
zero temperature ("quantum critical point") and asso- 
ciated effects due to quantum fluctuations at very low 
temperatures. Reference [l| gives an introductory review 
of the subject. 

These developments have generated a variety of dif- 
ficult theoretical questions; among them is the issue of 
how to treat the regime of strongly-interacting quan- 
tum fluctuations. In a recent paper Q, we developed 
an extension of the quasiparticle concept of Fermi liquid 
theory to the non Fermi-liquid regime near a quantum 
critical point (QCP). In essence, the theory goes beyond 
the Gaussian regime of critical fluctuations by introduc- 
ing interactions among the quantum fluctuations into the 
correlation function of the fluctuations. Central to the 
analysis is the concept of critical quasiparticlcs, which is 
based on the recognition that the single-particle spectral 
function can display a quasiparticle peak at non-zero ex- 
citation energy or temperature. This is expressed as a 
non-zero quasiparticle weight Z(w) for \u>\ not too small, 
although, as in a non Fermi liquid, at the Fermi surface 
Z(lo = 0) = 0. 

We realized the critical quasiparticle theory for the 
case of an antiferromagnctic quantum critical point and 
applied it to several quantities, principally resistivity and 
specific heat, for successful comparison to experimen- 
tal results on the heavy-fermion metal YbRh 2 Si2 (YRS), 
thereby showing that the theory, which describes a phys- 
ically transparent scenario, is capable of accounting for 
experimental results on a quantum critical metal. 

The implementation of the theory for an antiferro- 
magnctic (AFM) quantum critical point for a heavy- 
fermion compound is based on the recognition that be- 
low a "lattice Kondo temperature" Tkl, hybridization 
between conduction (s,p, d) electrons and local magnetic 
moments (/ orbitals) produces a heavy-electron liquid 
with an associated mass enhancement due to the orig- 
inally localized character of the / electrons. However, 



the / electrons are also responsible for the antiferromag- 
netism in a region of the phase diagram of the mate- 
rial. Near the AFM critical point, critical spin fluctu- 
ations are enhanced and interact with the heavy quasi- 
particles. This produces further mass enhancement and 
occurs in two stages. Above a certain temperature T x 
but below Tkl, the critical fluctuations are Gaussian (i.e 
non-interacting); those of two-dimensional AFM charac- 
ter (or three dimensional ferromagnetic character) pro- 
duce a logarithmic mass enhancement. A review of this 
physics can be found in reference Q. Below T x , how- 
ever, the critical fluctuations begin to interact with each 
other and consequently the effects on the quasiparticlcs 
change, e.g. the effective mass, or Z-factor, acquires a 
singular power-law frequency dependence. This is the re- 
gion described by the critical quasiparticle theory. Figure 
[1] illustrates these regions in the phase diagram for the 
heavy-fermion metal YRS. In YRS, an AFM quantum 
critical point is accessed by tuning a magnetic field H , 
but the theory is appropriate whatever the nature of the 
ordered phase and whatever the tuning parameter. 

Figure [T] shows the phase diagram for YbRli2Si2 near 
the quantum critical point. Five different regions may 
be seen: (1) The antifcrromagnetically ordered region 
at temperatures below the Neel temperature Tjy(iJ); (2) 
The Landau Fermi-liquid regime at T < Tfl{H); (3) 
A high temperature local- moment regime at T > Tkl, 
where Tkl is the characteristic temperature of the lat- 
tice Kondo effect, below which coherent heavy quasiparti- 
cles form by hybridization of /-electrons and conduction 
electrons; (4) A regime of quantum critical quasi- two- 
dimensional antifcrromagnetic Gaussian fluctuations at 
Tx ^$ T < Tkl, characterized by moderate non- Fermi 
liquid behavior; (5) The true critical regime, "QC", for 
T <T X governed by three-dimensional antifcrromagnetic 
fluctuations interacting strongly with the heavy quasi- 
particles. 

The basic input for the theory is a phenomenological 
form for the critical antiferromagnetic (AFM) spin fluc- 
tuation correlator at low temperature below T x , which 
reflects that it is generated from quasiparticles having a 
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FIG. 1. Phase diagram for YbRh2Si2 in the neighborhood of 
the critical magnetic field H c ~ 0.06 T: The dashed lines rep- 
resent crossovers. As the temperature is lowered below Tkl, 
the lattice Kondo effect and heavy quasiparticles develop and 
weakly-interacting ("Gaussian") 2d AFM fluctuations associ- 
ated with the quantum critical point cause non-Fermi liquid 
behavior. The T x line represents the crossover to strongly- 
interacting 3d fluctuations, which are dominant within the 
cone of quantum criticality (QC), bounded on the right by 
the crossover at Tfl{H) into the heavy Fermi- liquid state 
and on the left by the curved dashed line within the anti- 
ferromagnetic ordered state that sets in below Tn(H). The 
T*(H) crossover (dot-dashed line) is discussed in the text, 
below Eq. ([3j). 



non-zero frequency-dependent Z-factor: 



(N /Z)(co/v F Q) 



[r (H) + Z(q - Q) 2 £ 2 ] 2 + {lo/v f QY ' 



(1) 



where N is the bare density of states at the Fermi sur- 
face, v* F = (m,b/m*)vF is the renormalized quasiparticle 
Fermi velocity, £o — k F l is the microscopic AFM cor- 
relation length and r (H) = 1 + F(Q, H) oc H/H c - 1. 
Here F(Q, H) is a dimensionless generalized Landau pa- 
rameter, which — > —1 at the critical point. For conve- 
nience, we denote the underlying tuning parameter by 
H, which for YRS is the magnetic field. The interaction 
between quasiparticles and spin fluctuations / is related 
to to F(Q, H) by I = ZF{Q, H)/N . The "bare" quanti- 
ties nib, No, vf here are to be understood as those quan- 
tities already renormalized at Tkl by hybridization and 
the lattice Kondo effect, as discussed above. The addi- 
tional (frequency-dependent) mass enhancement caused 
by interaction with critical spin fluctuations is denoted 
as m* /nib = 1/Z. We do not consider questions rele- 
vant to the crossover Q from the local moment phase at 



T > Tkl to the heavy-fermion phase at T < T K l- The 
experimentally determined values in YRS are Tkl ~ 25 
K and T x « 0.3 K. 

The justification for the form of Eq. (JlJ is discussed 
in detail in reference We neglect the momentum de- 
pendence of the Z-factor by assuming that the critical 
behavior at the hot spots is spread over the whole Fermi 
surface by impurity scattering or other interaction effects. 
The Z-factor is determined by the quasiparticle self en- 
ergy, which in turn is determined by the interaction with 
the spin fluctuations. This leads [2J to a self-consistency 
relation in the form of a differential equation for Z(oj): 



= 1 + A— (Z 

duj 
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(2) 



This equation has two different physically meaningful 
solutions. The first is a weak coupling solution, valid 
provided the second term on the r.h.s. of Eq. @ is 
<C I in the energy range considered. This corresponds 
to the conventional spin-density wave scenario as dis- 
cussed in the works of Hertz @, Millis @ and Moriya 0. 
There exists, however, a second solution in the strong- 
coupling domain; it is accessible provided the initial Z _1 
at the scale when the three-dimensional antiferromag- 
netic fluctuation regime is entered (in the case of YRS at 
T « 0.3K) is sufficiently large such that the second term 
on the r.h.s. of Eq. (J2J dominates. We conjecture that 
in YRS, two-dimensional antiferromagnetic (assisted by 
three-dimensional ferromagnetic) spin fluctuations above 
T w 0.3K can provide the necessary growth of the effec- 
tive mass everywhere on the Fermi surface, even in the 
absence of impurity scattering. (In the case of quasi-2d 
antiferromagnetic fluctuations in a 3d metal this has been 
noticed first in reference 0]). Then the strong-coupling 
solution is 



3- i/2 (^or 3/2 KVQ) l/4 . 



(3) 



At non-zero temperature, Z is approximately given by 
Z(cj,T)oc[o; 2 + (7rT) 2 ] 1 /8. 

In the following, we describe the phase diagram ob- 
tained within our theory as consisting of the ordered 
phase bounded by the Neel temperature T N (H), the crit- 
ical regime ( "critical cone" ) at T x > T > Tp^ (H) and the 
Landau Fermi liquid regime at T < Tpp(H) (sec figure 
[TJ. In addition to these phase boundaries and crossover 
lines, Hall effect and other measurements in YRS 0-11 
show anomalies at a temperature T* (H) in the critical re- 
gion. The origin of the T*{H) line is not clear at present. 
Recent experiments on YRS under pressure fl2| or YRS 
doped with Co or Ir [l3| have shown that while the crit- 
ical field for the magnetic transition is found to shift, 
the T* line stays unchanged. This suggests that the T* 
feature is only weakly or not at all tied to the AFM tran- 
sition. The present theory ascribes the critical behavior 
to interaction of renormalized quasiparticles with renor- 
malized antiferromagnetic spin fluctuations, rather than 
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with fluctuations associated with T*(H), as has been sug- 
gested elsewhere [III EH Ell • 

In this paper, we elaborate the critical quasiparticlc 
theory further by deriving the scaling behavior of quan- 
tities near the QCP and we add the critical behavior 
of the Gruneisen coefficient, thermopower, thermal ex- 
pansion, magnetization, and susceptibility to the exper- 
imental consequences of the theory. Again, we compare 
these results of the theory to experimental information 
on YbRli2Si2, which has a magnetically tuned QCP. 

The deduction of the critical behavior of several ther- 
modynamic quantities from the scaling form of the free 
energy has been discussed by Zhu, et al [lij]. The corre- 
sponding exponents for these quantities are determined 
by the correlation length exponent v, the dynamical ex- 
ponent z and the dimensionality d. 



Scaling 

As discussed in reference , the critical regime is char- 
acterized by fluctuations of dimensionality d = 3. We can 
determine v and z from the critical quasiparticle theory 
as follows: It can be shown that when the frequency de- 
pendence of Z is a power less than 1/2, the frequency- 
dependent Z-factors appearing in Eq. ([!} have the argu- 
ment u>. Then from Eq. (Q]) and the result Z(lu) oc uj a , 
where we found a = 1/4, we can determine the dynami- 
cal exponent z = 4 since at the critical point 



Imx^a;)- 1 ~ Z 2 (io)q 4 + u 2 /Z 2 (uj). 



(4) 



Here and in the following, q denotes the deviation from 
the ordering wave vector Q. The scaling exponent of 
the spatial correlation length £ is found as v = 1/3 from 
r (H) ~ Zq 2 ~ {q z )^ 4 q 2 ~ q 3 so that £_ ~ r~ 1/3 . At 
T = we thus find £ oc (if — i7 c ) _1//3 .The characteristic 
correlation time is given by t c ~ £ 2 ~ (H — H c )~ 4 / 3 . 

When H = H c {i.e. r = 0), the approach to the 
T = critical point is described by a temperature de- 
pendence of the control parameter, which can be written 
r{H Cl T) ~ T x l zv . The condition r(H c , T) « r(H, T = 0) 
or T = T FL (H) oc \H-H C \ ZV defines the boundary of the 
"critical cone" of quantum criticality. Then, within the 
quantum critical regime, the temperature dependence of 
the correlation length is found as £ oc T~ x / Z . It may 
be shown that the above critical temperature dependent 
correction to r is caused by the interaction of spin fluc- 
tuations treated at the Hartree level [2lJ, giving rise to 
r{H Cl T)~T/Z{T). 

From these results (v = 1/3, z = 4), we may write 
Imx(q, oj), or equivalently the structure factor S(q,u>), 
in the scaling form 



S(<?,w) = £ 4 $K 4 , 9 0, *(x,j/)oc 



„l/2 



Note that at the critical wave vector (at q = 0), uj/T 
scaling holds. 

The scaling form for the free energy density may be 
deduced by dimensional analysis 0] : 



f(H,T)=r id+z) ^f(r o e /u ,T^) 



T~ 



'0 



(6a) 
(6b) 



In our case, d — 3, v = 1/3, z = 4. 

We may generalize these results to arbitrary dimension 
d. First, we note (from the derivation of the quasiparticle 
weight factor Z presented in reference 1) that d = 4 is the 
upper critical dimension, so that the spin fluctuations do 
not destroy the Fermi liquid state for d > 4. At d < 2, we 
find that the critical quasiparticles are no longer well de- 
fined. We therefore consider here only dimensions within 
the interval 2 < d < 4. Using the exponent a to describe 
the frequency power law of the quasiparticle weight Z, 
we find, from Eq. (|3J), that the dynamical critical expo- 
nent is given by z = 2/(1 — 2a) and the correlation length 
exponent is v = 1/(2 + za) 



Specific heat 

The critical part of the specific heat at the criti- 
cal tuning parameter r = is found from Eq. (|6aj) : 
C c oc T d l z = T 3 / 4 as we found before 0. We can obtain 
the behavior of the specific heat near T = as a func- 
tion of r from Eq. (|6bl) as follows [l(|: Near x = 0, ^Sf(x) 
behaves as $(x ->■ 0) = *(0) + bx**. Therefore 

C(r,T) oc r^+^'^^i - 1)T" -1 

= ^(/i-l)T" -1 r- 1/3( ^ _7) . (7) 

This behavior is valid in the region of the T, H phase 
diagram defined by T jr zv < 1, which is T < Tfl(H) oc 
\H-H C \ Z ». 

As we showed in reference Q , the theoretical behavior 
of 7 = C c /T at criticality found above agrees very well 
with the data. For the low T-dependence of 7 on r = 
\H — H c \, we note that in YRS on either side of the QCP, 
Fermi liquid behavior obtains ljj so that /1 = 2. Then, 
from Eq. [7J 7 = C/T varies as \H - F c |" 1/3 . This is 
precisely the behavior observed in experiment on YRS 

ilGl- 



(l + a; 1 /4 2/ 2)2 +;z .3/2- 
(5) 



Resistivity 

The scaling analysis may be applied to determine the 
temperature and magnetic field power-law dependences 
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of the resistivity, which (assuming impurity and umklapp 
scattering) is given by: 

p(T,H) - p(0) <xT(T,H)/Z(T,H), 

where T is the quasiparticle relaxation rate, which we 
found [2] at d = 3 and H = H c , to be proportional to 
ufil 2 jZ 2 (uj). Since T has dimensions of energy, its scaling 
form is r(T, H) = £- z $ r (T£ z ) as in Eq. With 
Z(oj) oc oj a oc £,~ za , we have 

p(T, H) - P (0) oc e (Q_1)z $r(T^ 2 ) 

To insure that p is finite at H = H c and T > 0, we 
require <f>r{z) oc x x ~ a as x —¥ oo, so that 



p(T,H c ) -p(0) ocT 1 



(8) 



As we discussed above for the specific heat, away from 
criticality, at low temperature outside either side of the 
critical cone defined by T = Tfl, Fermi liquid behavior, 
p{T) - p(0) oc T 2 obtains, so that $ r (x) oc x . Then 



p(T,H) -p(0) oc£ (q+1)z T 2 . 



(9) 



Using the values a = 1/4, z = 4, which we found in the 
self-consistent critical quasiparticle theory Q, we find 
that Eqs. (6,7) become 

p(T,H c )-p(0) ocT 3 / 4 (10a) 
p(T, J3") - p(0) oc £ 5 T 2 oc |i? - H c \- 5/3 T 2 . (10b) 



The available data agree very well with Eq. (|10a[) as found 
in reference 0] by direct calculation. At the time of this 
writing, for the field dependence of the coefficient of T 2 
in Eq. (|10b[) . there is insufficient data close enough to H c 
to allow a quantitative comparison. 



Thermopower 

At very low temperatures, when impurity scattering 
dominates, the thermopower S may be expressed in terms 
of the electrical conductivity <7i m p(p), depending on the 
chemical potential p, by the Mott formula 



s _ if_ T d In a imp 



e 



dp 



(11) 



Since the mean free path £ is not renormalized, 
Oimp{p) oc kp£, where &f is the Fermi wave vector, is 
independent of the effective mass. The derivative with 
respect to p may be expressed by the one with respect to 
the Fermi wave vector as d/dp = (\/v* F ){d / dkp) and 
ddimp/ dkp is independent of m*. Therefore, we find 
S oc Tm*{T). In the Gaussian fluctuation reg ime, we 
get S oc Tln(T), which has been observed [18(. In the 
critical regime proper, at the critical field H = 0.064 T 



we predict S oc T 3 / 4 . The comparison of our result to 
the experiment reported in reference 18( is shown in fig- 
ure [21 Unfortunately, the data was taken at H = 0.06 T, 
below the critical field of 0.064 T. At the measurement 
field, the AFM Neel temperature T N is about 20 mK [lH, 
so we do not expect good agreement at the very lowest 
temperatures, as seen in the figure. 




T(K) 

FIG. 2. Thermopower: Comparison of theory and data of 
reference [l^| near the critical magnetic field and below the 
critical temperature for quantum critical scaling. 



Magnetization and susceptibility 

The temperature dependence of the magnetization 
may be found from the derivative of the free energy in 
Eq. (paj) with respect to field: 



df_ 
dH 



1 fdf 



T H c \dro /T 



(12) 



Assuming M is regular at ro = so that ^'(0) = const, 
we get 

.d+z-Uu 



M(T) - M(0) oc — T" 



-T. 



(13) 



This result agrees fairly well with the experimental data 
at low T (below 400 mK) as reported in reference . A 
very good fit is obtained if the critical T-dependence of 
Eq. (fTB"]) is augmented by the sub-leading T 2 contribution 
due to Fermi-liquid effects. This is shown on figure [3j 
where the theoretical behavior at H = 0.06 T is given by 

M(T) - M(0) = -0.68T + 0.63T 2 . (14) 

For the field dependence at small T, we use Eq. (|6b[) : 



Of 



1 



M crit = -(^J =-—v(d + z)r Q 



v{d+z)-l 



T=0 



*(0), 
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which for d — 3, z = 4, v = 1/3, ro = H/H c — 1 is 



M(H) - M(H C ) 



7 1 / H 



3 H r \H. 



4 / 3 „ 
1 I tf(0). (15) 



Therefore, the critical field dependence of the suscepti- 
bility is 

Xo I it=g W <x-(H-H c ) 1 ^. (16) 

The available data [itj indicate a cusp-like behavior of 
x(H) when approaching H c , but do not extend close 
enough to allow for a quantitative comparison. 



M = 0.39 - 0.68T + 0.63T 2 
(at H = 0.06 T) 




0.15 0.2 
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FIG. 3. Mag netization: Comparison of theory and data of 
reference [l!| near the critical magnetic field. A T 2 Fermi- 
liquid contribution has been added to the critical behavior. 




FIG. 4. Magnetic susceptibility: Comparison of theory, Eq. 
(|17|l and data of reference [52] near the critical magnetic field 
and below the critical temperature for quantum critical scal- 
ing. Similarly to Eq. (fT4")> and figure [3] we have added a T 2 
Fermi-liquid contribution. 

The critical T-dcpcndence of the susceptibility is found 



from Eq. (|T2J) as 

dM 
~dH 



1 



\ 



dM i±* 
cx T 



-2/. 



H c dr 



T 1 / 4 , (17) 



where we assumed ^"(0) = const. We can compare this 
result to the experiment reported by Gegenwart et al [22J , 
where %ac was measured below 0.2 K at a field of 0.05 T. 
Since Tn is about 30 mK at this field, we do not include 
data below 40 mK in the comparison shown on figure |4l 



Griineisen ratio 

In the case of a magnetic field tuned QCP, as in YRS, 
the magnetic Griineisen ratio fltij is studied: 



r 



JVf 



{dM/dT) H {dS/dB)i 



(18) 



The critical part of the Griineisen ratio has been dis- 
cussed in reference fl(| and measured in YRS by Tokiwa, 
et al [ijj . As in reference [l6| , the behavior of Tm may be 
determined from the scaling analysis. We find, for z = 4, 
v = 1/3 and n = 2, that [dS/dH] ra=0 cx T°. Then 

-3/4 



r M (r = 0,T)«T- 
r M (r,T^0) = 



-1/3 



H - H r 



(19a) 
(19b) 



The factor — 1/3 = v(d — z) in Eq. (QjjbJ is universal as 
the scaling functions in the ratio defining V m cancel. 

These results as well as those for the thermodynamic 
Griineisen parameter can also be obtained [2l| by direct 
calculation from derivatives of the self-energy expressions 
given in reference Q. 

Just as for the specific heat and the resistivity 0], 
there is excellent agreement between the critical quasi- 
particle theory and the experiment [l9| for the magnetic 
Griineisen ratio. The factor v(d — z) = —1/3 in Eq. (|19b[) 
is measured as —0.30 ± 0.01 and the comparison of our 
result of Eq. (|19ap of the T~ 3/4 dependence at H = H c 
to the experiment is shown on figure [5] 

We may evaluate the consistency of our results by using 
the theory fits for M (T) , Eq. ([T4"|) and for the specific heat 
shown on figure 1 of reference @ (C = 0.034 + 0.454T 3 / 4 ) 
in Eq. (|T8)) . There are no further adjustable parameters 
and the agreement is shown in the inset to figure [S] 

The volume thermal expansion coefficient P(T) enters 
the conventional thermodynamic Griineisen ratio Tt = 
(3(r ,T)/C(r , T). The critical behavior here is tuned by 
pressure P, so that r oc (P c — P)/P c . The critical part 
of /3 is determined by derivatives of the free energy as 



Pc{r ,T) 



1 d 2 f 
Pv dTdr Q 



ocT°, 



where v is the specific volume and we have used Eq. (|6aj) . 
To this we may add the conventional Fermi-liquid con- 
tribution cx T so that j3 = a + bT in the critical region. 
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FIG. 5. Magnetic Griineisen ratio: Comparison of theory, 
Eq. (|19a[) and data of reference Q1J at the critical magnetic 
field and below the critical temperature for quantum critical 
scaling. The inset shows the comparison of the data and Eq. 
(I8jl when the previously fit M(T) and C(T) are used. 



Conclusion 

The critical quasiparticlc theory describes a renor- 
malized Gaussian picture of critical fluctuations, where 
the fluctuations arc assumed to be non-interacting be- 
yond effects that renormalize the quantities entering the 
fluctuation spectrum. The renormalizations of these pa- 
rameters, the staggered static susceptibility and the Lan- 
dau damping term, are, however, of a strong-coupling 
nature. It remains to be shown that fluctuation interac- 
tions are irrelevant within a more complete theory. We 
have used a scaling analysis of the theory to derive the 
critical behavior of a number of experimentally observed 
quantities in the heavy fermion compound YbRli2Si2. 
The good agreement with all these observed quantities 
indicates that the critical quasiparticle theory captures 
essential features of the critical behavior in this com- 
pound. We note that YbRh 2 Si2 is especially suitable 
for application of the theory since it has the wide region 
of two-dimensional (2d) critical behavior that is neces- 
sary to ultimately access the strong coupling regime of 
three-dimensional fluctuations at lower temperatures, as 
discussed in reference , where the critical quasiparticle 
theory was introduced. 

Besides YRS, other candidate systems have been stud- 
ied with different degrees of detail. For CeCue ^Au^ 
it has been found (for a review see reference [3() that 
some of the low-temperature properties like the specific 
heat and the resistivity indicate quasi-2d antiferromag- 
netic fluctuations in the Gaussian regime. Neutron scat- 
tering studies at somewhat higher energies have been in- 
terpreted as showing critical scaling behavior incompat- 
ible with the Gaussian theory. It is conceivable that the 



behavior is different depending on the energy scale. Fur- 
ther candidate systems like CeCoIns also appear to show 
quasi-2d antifcrromagnetic Gaussian fluctuations (for an 
overview, see reference Q). In all of these systems we ex- 
pect a crossover to 3d fluctuations at lower temperatures, 
in which case our theory should apply. 
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